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ABSTRACT 

Essentially all stars form in giant molecular clouds (GMCs). However, inside GMCs, most of the gas does 
not participate in star formation; rather, denser gas accumulates in clumps in the GMC, with the bulk of the 
stars in a given GMC forming in a few of the most massive clumps. In the Milky Way, these clumps have 
masses M Q \ <5x 10~ 2 of the GMC, radii r c i ~ 1 pc, and free-fall times r C ] ~ 2 x 10 5 yr. We show that clumps 
inside giant molecular clouds should accrete at a modified Bondi accretion rate, which depends on clump mass 

as M c i ~ M 5 ( A . This rate is initially rather slow, usually slower than the initial star formation rate inside the 
clump (we adopt the common assumption that inside the clump, = egM c ijT c i, with eg ~ 0.017). However, 
after ~ 2 GMC free-fall times tgmc> me clump accretion rate accelerates rapidly; formally, the clump can 
accrete the entire GMC in ~ 3tgmc- At the same time, the star formation rate accelerates, tracking the Bondi 
accretion rate. If the GMC is disrupted by feedback from the largest clump, half the stars in that clump form in 
the final tgmc before the GMC is disrupted. The theory predicts that the distribution of effective star formation 
rates, measured per GMC free-fall time, is broad, ranging from ~ 0.001 up to 0.1 or larger and that the mass 
spectrum of star clusters is flatter than that of clumps, consistent with observations. 

Subject headings: galaxiesTSM — galaxies: star clusters: general — ISMxlouds — stars: formation — turbulence 



1. INTRODUCTION 

The star formation rate (SFR) is a fundamental parameter of 
disk galaxies. It is well characterize d on galactic disk scale s 
by the Kennicutt-Schmidt relations dKennicutlll 19891 119981) . 
These relations come in two forms. The first form is that of a 
correlation between the surface density of star formation 
(in solar masses per year or grams per second) and that of gas, 

^gas> 

£*=A£« S , (1) 

with aw 1 .4. The second form relates to the surface den- 
sity of gas via the disk dynamical time £1 = v c /Rd, 



(2) 



where Rj is the disk radius and v c is the circular velocity 
of the galaxy. The dimensionless parameter has an observa- 
tionally determined value r\ « 0.017. More recent work has 
refined these relat i ons, particularly a t low surface densities 
dLeroy et al J 120081: iBigiel et al J 120081) . but in galaxies where 
the bulk of the gas is molecular, Equations (0 and remain 
valid. 

There are numerous theoretical explanations of these re- 
lations, relying on ver y different physics , including sup- 
port by magnetic fields dMouschoviasll 1 9761 Shu 1983 b, sup- 
pression of collapse by sup ersonic turbulence (Padoan Hl995l: 
Krumholz & McKee 2005), energy feedback from stars and 
supernovae, and m omentum feedback from massive stars 
dMurrav et al. 2010). Which of these explanations is correct, 
if any, is currently still under debate. 

Other work has focused on relating the rate of star for- 
mation to gas surface or volume density on sub-disk scales, 
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ranging from ~ 1 kpc ( Schruba et al. 20 To|) in nearby galax- 
ies, to ~ lOOpc in the Milky Way dMoonev & So lomon 
1988[lMead et al.ll990tlEvansll991tlLada et al.l2010tlMurrayl 
201 ll) . These studies find relations similar in form to Equation 
10. The mean value of the coefficient, averaged over many ar- 
eas or star forming regions, is consistent with the global value 
rj ss 0.017. However, the dispersion of r\ appears to vary with 
the scale on which the star formation is probed: measured 
values ranges from r\ < 10~ 3 to r] ~ 0.5 - ranging almost four 
orders in magnitude. 

It is well established that not all the gas in a galaxy par- 
ticipates in star for mation; star formatio n takes place only in 
molecular gas, ("e.g. lSchruba et al.ll201 ll) . Further, even inside 
a giant molecular cloud (GMC), not all the gas participates; 
rather, stars form primarily in high de nsity gas, often in th e 
form of "clumps" and filaments (e.g. Mol inari et al.1 120 lOh . 
Another way to see this is that star formation rates scale non- 
linearly with CO luminosity, which traces rather low density 
gas, while t hey scale linearly with HCN, which traces high 
density gas dGao & Solomon! [20041) . In local spirals like the 
Milky Way, the fraction of HCN gas relative to mole cular (CO 
emitti ng) gas is roughly several to ten percent (e.g IWu et al.l 
l2OT0h . 

In this paper we present a simple theory for the rate of 
star formation in GMCs, in which gravity plays the domi- 
nant role. We assume that the gas in a GMC is turbulently 
supported, i.e., the GMC is in rough virial equilibrium. The 
turbulence seeds the cloud with parsec scale clumps hav- 
ing masses 6 = M c1 /Mgmc ~ 10" 3 ; IWu et al.l d2010b refer to 
these objects as "massive dense clumps". These clumps then 
grow by accretion. While feedback from one or two of the 
most massive clumps eventually overcomes gravity and dis- 
rupts the GMC hosting the star fo rmation (see for instance, 
lMurravetaill2010t lMurravll20rTl) . we ignore the effects of 
feedback. We argue that the feedback affects the accretion 
only in the last stages of clump growth; both the clump accre- 
tion and star formation rate accelerate rapidly, and hence so 
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do the effects of feedback. 

By analogy with Equation (|2), we define the effective GMC 
star formation efficiency per free-fall time as 



M., 



»7ff,G = r GMC 



M G MC 



(3) 



where tgmc is the free-fall time of the GMC. Our gravity- 
dominated theory of star formation in GMCs predicts that r}n q 
is a strong function of time; it is small immediately after the 
GMC is assembled and for a substantial fraction of a GMC 
free-fall time, but then increases rapidly; we show that this 
explains the large dispersion of r\ on small scales referred to 
above. 

This paper is organized as follows. In ^2] we describe ac- 
cretion onto massive clumps in GMCs; the mass accretion rate 
accelerates as a clump grows. We show that the star forma- 
tion rate in the clump tends to track the mass accretion rate. It 
follows that rjs t o is time dependent, even if the star formation 
rate per free fall time (e^, defined in Equation (O below) is 
constant for the parsec scale clumps of molecular gas which 
form individual star clusters. In ^3]we compare our results to 
observations of star formation in the Milky Way and nearby 
galaxies; the theory predicts a large spread in the apparent star 
formation rate in GMCs, and a flatter mass distribution of star 
clusters compared to star forming clumps. In $4]we compare 
our results with recent numerical simulations, and we briefly 
discuss the effects of vorticity, arguing that they do not limit 
the rate of accretion on to dense clumps. Finally, we close 
with our conclusions. 

2. BONDI ACCRETION IN GIANT MOLECULAR CLOUDS 

Simulations of supersonic turbulence suggest that it pre- 
vents gas from rapidly fra gmenting and turning into stars on 
the local dyna mical time (Padoan 1995': iKfessen et ai]l2000t 
iLi et al.ir2004l) . Based on these and similar results, numerous 
authors have suggested that, over the local dynamical or free 
fall time, only a s mall fraction of the gas in a given regio n 
is turned into stars (Padoan 1995; Krumholz & McKee 2005). 
In particular, if the mean density in a region is p, the local 
free-fall time is 



Tff ; 



/ 3tt 
32Gp' 



Then the star formation rate is given by 



dM* 
dt 



■ £ff- 



Mi 

Tff 



(4) 



(5) 



with £ff w 0.017. While we are agnostic about the applicabil- 
ity of this prescription, we will adopt it here. 

We now show that the local star formation rate given by 
Equation (j5]l does not limit the overall rate of star formation in 
GMCs very significantly. This results from two facts; first, for 
dense enough clumps, the star formation timescale is shorter 
than eff times the free-fall time of the GMC, and second, star 
forming clumps control the dynamics of the gas in their vicin- 
ity through their gravity, as we now show. 

GMCs are observed to harbor both massive gas clumps and 
massive star clusters. For example, in the Milky Way, the 
most massiv e GMCs have Mqmc ~ 3 x 10 6 M^ and R GM c ~ 
lOOpc, e.g.. iBronfman et all dl988h : iGrabelskv etaLl (1988). 
These massive GMCs contain dense clumps with M c \ ss 
10 4 M o dSimon et al.1 1200U I Wu et all 120101) and star clus - 
ters with similar or larger masses jMurrav & Rahmanll2010l) . 



We show that the masses of these gas clumps initially grow 
slowly, but the growth rate then accelerates rapidly, so that, 
left unchecked, the largest clump would consume the GMC 
in a few Tff. 

Consider a giant molecular cloud of radius Rcmc and mass 
Mqmc> which we will assume is near virial equilibrium. The 
velocity on the scale R CMC is v GM c = \/GM CM c/Rgmc, while 
the turbulent velocity inside the cloud is given by 



vr(r) = vqmc 



Rg 



MC 



(6) 



where p ~ 1/2, by one of Larson's laws dLarsonl I198U 
ISolomon et al ll987l) . In the most massive Milky Way GMCs, 
^y(^gmc) = v'gmc ~ 6kms _1 ; in more rapidly star forming 
galaxies such as the z = 2 galaxy BX 482, vgmc ~ 50 km s~ . 

Now consider a clump embedded in the GMC with a total 
mass M c \ (including a stellar mass M„) and radius r c \. The 
gravity of the cluster co ntrols the flo w of the surrounding gas 
out to the Bondi radius ( Bondi| [T952l) 



GM cl 

V'tOb) 



Solving for r B , 



T7 "GMC 

VMgmc/ 



M 1/2 *g: 



MC, 



(7) 



(8) 



where we have defined p = M c \ / Mqmc ■ The turbulent velocity 
at the Bondi radius is 



VrOB) = M 1/4 V GMC- 

The associated Bondi accretion rate is 

M B = 4nrlv T (r B )p. 



(9) 



(10) 



The density p is a function of position inside the GMC. As- 
suming spherical symmetry, the literature contains suggested 
functions of the form 



p(r) = Cp(^ 
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where C = (3 - f3)/3, and p = 3Mgmc /(4t^gmc) * s trie mean 
density of the GMC. We will consider the case of /3 = (con- 
stant density), 1, and 2 (isothermal distribution). 

It is likely that clumps form where gas flows shock and sub- 
sequently cool. While such shocks can occur anywhere in 
the interior of the GMC, there is some observational evidence 
that star clusters may be biased toward the centers of GMCs, 
e.g., HII regions are c entrally concentrated in their host GMCs 
dScoviile et al.| [l987). We will simply assume that our clump 
forms at a distance R c i from the center of the GMC. 

2.1. Constant density GMCs 

For a constant density cloud (/3 = 0), the position of the 
clump does not affect the accretion rate (so long as the Bondi 
radius of the clump does not exceed the distance to the surface 
of the GMC) so we will focus on this case first. The accretion 
rate is then 



Mi 



Bondi " 



3tt 
2V2 



5/4 



MqmC 
T GMC 



(12) 



where tgmc is the free-fall time of the GMC. 
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Milky Way GMCs often have m any such clumps o ver a 
ranges of masses (see for instance ISimon etail [2001 ). We 

5/4 

focus on the largest clump because the M J dependence on 
the growth rate in Equation ( TTZb implies that the largest clump 
also grows the fastest. We assume that clumps are born with 
an initial mass Mq = <5effMcMC, and that the sum of all the 
gravitationally bound clumps produced over a free-fall time 
is ~ EffMcMC, i.e., that the sum over clumps Y,5 = 1. 

The time evolution of the mass of an individual accreting 
clump is given by 



M cl (r) = M 



3tt 



Mq 

8\/2tgmc V m gmc 



1/4' 



-4 



(13) 



The denominator in Equation ( fT3l vanishes at a finite time; 
before that happens, the GMC is consumed, after a time 



ffinal = 



8\/2 
3^T 



/ M CMC \ 
V M {) ) 



1/4 



1 



'3.5 



_s_ 

02 



-1/4 



eg 
0.02 



-1/4 



TGMC- 



(14) 



(15) 



We interpret ffi na i as a hard upper limit on the lifetimes of 
GMCs; it is the time it takes the clump to completely ac- 
crete the GMC. While this calculation ignores both angular 
momentum considerations and feedback processes, we argue 
below that the GMC will be disrupted by stellar feedback be- 
fore ffinai is reached and that angular momentum does not play 
a significant role. 

2.1.1. Star Formation in Bondi Clumps 

The discussion in the previous section follows the evolution 
of a clump as it accretes gas from its parent GMC over several 
tgmc- Now we refine our picture, and discuss the evolution of 
stellar mass in the clump. The clump mass is divided into 
stars and gas; we denote the stellar mass as M*, and we define 
the gas fraction 

t M cl -M* 

f e = ^T- (16) 
We use the star formation law from Equation ©, 



dM* 
dt 



■■ e ff - 



fgMg 

Tel 



(17) 



where t c \ is the free-fall time of the cluster. The free-fall time 
of the cluster is given by 



T 1 = 



3tt 



32Gp c 



= tgmci 



I(Rci/Rgmc) 3 



(18) 



It is useful to re-scale all the cluster masses (total, stellar, 
and gas) by the GMC mass, so in addition to \i = M c i/Mgmc, 
we have ^ g = f g [i and [i* = M*/Mgmc- The evolution equa- 
tions of the gas and the stars are then given by 



Ms 

M* = £ ff — > 

Tel 

^g = /i-e ff — . 

T C 1 

Combining Equations ( TToT ) and d20l > we have 
— - i l -Jg)-- e «— 



dt 



/' 



Tel 



(19) 
(20) 

(21) 



We can look for a fixed point of the last equation by setting 
the left-hand side of this equation to zero: 



f g 



1 



fixed ' 



i +*(/*)' 



where 



, , 2V2 fR CMC \ 



3/2 



3ir \ r c i(/x) / 



1/4 



(22) 



(23) 



We have called fg,& xe d a fixed point, although it is only fixed 
if is a constant. In fact, jl/ji ~ /i 1//4 is an increasing func- 
tion of time, so Equation (|221 does not give the instantaneous 
value of f g . As the cluster grows /i will increase, resulting in 
a change in f g *. The actual value of f g will be slightly differ- 
ent than / g «, with the magnitude of the off-set controlled by 
the ratio of the star formation time (from Equation ©) to the 
clump accretion time (from Equation (fTUll). 

It is easy to show that this psuedo-fixed point is an attrac- 
tor, in the sense that if f g is initially smaller than / g * the gas 
fraction will increase, approaching / g * asymptotically, while 
f g will decrease if it is larger than f gtf . 

2. 1.2. The mass radius relation for clumps 

In order to integrate Equations ( TTOb or (fT2l . ( fT9l l. and (l20l i. 
we need to know t c \, which requires knowledge of the mass- 
radius relation for massive clumps in GMCs. We have not 
been able to find much information in the literature, so as a 
proxy we use the mass-radius relation for young star clusters, 
under the assumption that the radius of such clusters reflects, 
to some extent, the radius of the parent gas clumps. 

The half light radii of star clusters (IWalcher et afl 120051: 
[Murray 2009; Scheepma ker et al.l2009l) and the effective radii 
of clumps (Si mon et al.l 12001) have characteristic values of 
1 -3pc for masses below 1O 6 M so we adopt a fiducial value 
of lpc. More massive star clusters have radii given by (see 
lMurravll2009h 



r c i(M d ) : 



if 



V !0 6 M f 



3/5 



pc. 



(24) 



These relations are clearly very rough approximations; bet- 
ter observational data on the mass-radius relation for massive 
clumps in GMCs would be of great value. 

Taking parameters from observed massive star clusters, we 
find that g(/i) ranges from a minimum value of g(fi) « 0.1 
for low mass (10 3 M Q ) clusters in 30- lOOpc, 10 5 - 1O 6 M 
GMCs in local galaxies up to a maximum value of g ps 50, for 
1O 6 M clusters in the R G mc = lkpc, 1O 8 M GMCs in high 
redshift star forming galaxies. The corresponding limits for 
the gas fraction are 0.02 < f g fi xe d ^ 0.9. 

2.1.3. Low mass star clusters in local galaxies 

The result of integrating Equation ( |2TT > for a cluster with 
Mci < 10 6 M Q in a 10 6 M Q , ^ G mc = 100pc GMC is illustrated 
in Figure Q] This shows the result of two numerical integra- 
tions, one with an initial gas fraction of one, and a second 
with an initial gas fraction of zero. The dashed line shows 
/g.fixedGu); me f act mat both numerical integrations converge 
to a value slightly above / s .fi xe d shows both that Equation (fJTJ 
has an attractor, and that / gl flxed is a good approximation for 
that attractor. Since r c i is taken to be fixed, g(p) is mono- 
tonically increasing, but it never exceeds unity, so f g {p) is 
monotonically decreasing and rather featureless. 
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FIG. 1 . — The gas fraction of an accreting clump plotted as a function of the 
ratio of clump mass to GMC mass, /i. The upper solid line shows the result 
of starting with a pure gas clump, while the lower solid line shows the result 
of starting with a purely stellar clump. The two lines rapidly converge as the 
clump mass grows. The dashed line sho ws t he approximate fixed point gas 
fraction f g * given by Equations )22t and 423K which applies in the limit of a 
very slowly growing clump fi < < 1 . 

Note that f g varies slowly with /i, only by a factor of 2 (af- 
ter the initial transient), while /i varies by nearly two orders 
of magnitude. This shows that the star formation rate traces 
the clump mass accretion rate. In other words, the star forma- 
tion rate is set by the clump accretion rate rather than by the 
properties of the local turbulence. 

2.1.4. Massive star clusters 

ForM c i > 1O 6 M , the clump radius grows with increasing 
mass. In this case, decreases with increasing fi, so that 
f g increases. Hence, the star formation rate lags behind the 
Bondi accretion rate. However, 1 — f g « gi/j) is never very 
small, as noted above, so the star cluster mass is always a 
substantial fraction of the clump mass - both increase rapidly 
as t —> ffi na i. Once again, the star formation rate is set by the 
clump mass accretion rate rather than by turbulence. 

In Figure |2j we plot the evolution of a Bondi clump accret- 
ing from a constant density GMC, where Mgmc = 1O 7 M and 
Rqmc = 100 pc. We show the total (solid line), stellar (dashed 
line) and gas (dotted line) mass fractions for the Bondi clump. 
For most of the lifetime of the clump, it remains small. How- 
ever it grows rapidly at late times, consuming the GMC after 
~ 4tgmc ~ 19Myr, in line with our expectations for ffi na i. We 
also track the stellar and gas mass separately and note the dif- 
ference in behavior for //* < 0. 1 (M* < 10 6 M Q ) and //* > 0. 1 
(M* > 10 6 M Q ). Initially the stellar mass tracks the clump 
mass very closely. Beyond M* > 10 6 M Q the stellar mass still 
tracks the clump mass, but not as closely as in the fixed r c \ 
regime. This is especially evident in the kink in ii g . 

2.2. Clump growth by Bondi Accretion: General Case 

We now consider the case of a clump located away from 
the center of a non-uniform GMC. Suppose the clump is a 
distance R c \ away from the centerQ If R c i ^> r B , then our pre- 
vious discussion holds - the clump accretes from a (locally) 

4 Note the difference between R c \ and . R c i defines the radial position of 
the clump in the GMC, while r c i defines the size of the clump. 



FIG. 2. — The total (solid line), stellar (dashed line), and gas (dotted line) 
mass fraction of a clump undergoing Bondi accretion, plotted as a function 
of t/TQMC- The mass and radius of the host GMC are Mgmc = 1O 7 M0 and 
Rqmc = lOOpc. The dynamic time of the GMC is tqmc = 4.5Myrs. We 
ignore the effects of stellar feedback. The clump experiences a long period 
of quiescience, which is followed by a period of rapid growth. The stellar 
mass initially tracks the clump mass very closely, but for M t > 10 6 Mq (at 
t/ T GMC ~ 3.3), the star formation rate grows less quickly than the accretion 
rate, so the stellar mass stops tracking the clump mass so closely. This is most 
easily seen by the change in slope of fJ. g (t). This change in behavior is due 
to the change in the clump mass-radius relation, which we have assumed is 
constant for star clusters with M* < 10 6 Mq, but has radius increasing with 
mass for M„ > 10 6 M Q . 



nearly uniform medium. If we relax this assumption, the den- 
sity depends on the position at the surface of the Bondi sphere. 
Recall that we assume a density profile of the GMC given by 

p(r) = po (r/RcMc) ■ 

We choose a second coordinate system centered on the 
clump, with the center of the GMC located on the z-axis a 
distance R c \ away. A point on the surface of the Bondi sphere, 
described by the second coordinate system (r B ,9, <fi) lies a dis- 
tance r(r B ,9, 4>) from the center of the GMC, where 

r(r B ,e,4>) = ^rl+Rl x + 2r B R cX cos 9. (25) 

We assume that the turbulent velocity is homogeneous and 
isotropic, unlike the density. Using this assumption, the mass 
accretion rate is: 

M= J rldcos6d(j)p(r(r B ,0,(j)))rlv T (r B ) 

(3-^GMC, 



5/4 



d cos 9 (R^ + Rl + 2R ci R B cos ( 



-13/2 



(26) 



where R B = r B /R CM c = y/fl, ^gmc = t g ^ c , and R ci = 
Rci/Rgmc- Performing the integral over cos 9, we find 



. (3-/3)S!gmc 3/4 

u= — a ' x 

2(2 -flRa P 



(« c 2 , + M+2^c1 V / m) 



1-/3/2 



(Rl + H-2R dy /ji) 



1-/3/2 



(27) 
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for j3 y 2 and 



(3-/3)ft GM c 3/4, / ^ c 2 1 + ^ + 2 ^c1 1 /M 
-jU ' log 1 



2fl r 



R 2 cl + H-2R cl ^jl 



(28) 



for /3 = 2. 

Equations (|27T i and d28l ) define the evolution of an accreting 
clump at a position /? c i. We note that Equation ( f27T > reduces to 



^ 3a p n GM cM 5/4 ^f 



1 + <9(/j 3/2 ) 



(29) 



for r B << /? c i, which we recognize as Equation (fT2l . where 
the background density is set by the position of the clump. A 
similar limiting form is found from Equation d28l . 

It is straightforward to integrate Equations (|27T i and (l28l . 
along with Equations ( [191 ) and ( [20) . The results are shown 
in Figures [2] (for p(r) = const., i.e., /3 = 0) and [3] for /3 = 1 
(left hand panel) and (3 = 2 (right hand panel). Here, we show 
the total (solid line), stellar (dashed line) and gas (dotted line) 
mass fractions for the various Bondi clumps. For /? ^ we 
present integrations at three example radii, R c i = 0.1 (black 
lines), 0.3 (blue lines), and 0.5 (red lines), to demonstrate that 
the clump mass accretion rate depends rather sensitively on 
the position of the clump in the GMC0 

Clumps in centrally concentrated GMCs exhibit rapid 
growth. This is shown by the R c \ = 0.1 curves (solid lines) 
of Figure [3] As we have already discussed, the high densities 
of these central regions lead to rapid growth, in which the en- 
tire GMC is accreted in a single (mean) GMC free-fall time 
tgmc- 

Clumps that are significantly off-center, i.e., R c i = 0.3 and 
0.5, behave more like the constant density GMC case. There 
is a period of slow growth, followed by a period of rapid 
growth. However, because these host GMCs are more cen- 
trally concentrated than a constant density GMC, the time it 
takes an off-center clump to accrete a substantial fraction of 
the host GMC mass is less than that seen in the constant den- 
sity GMC by a factor of a few, when measured in units tqmc- 

The division between "small" and "large" depends on 
the (initial) Bondi radius of the clump; since we are using 
/i = 0.003, the initial Bondi radius is 0.06, so R c i = 0. 1 is at the 
transition between rapidly accreting clumps and more slowly 
accreting clumps. 

3. COMPARISON TO OBSERVATIONS 

3.1. Bondi Clumps and the observed rate of star formation in 

GMCs 

These results suggest a robust story for star formation in 
GMCs, independent of the host GMC density profile. First, 
much of the star formation happens in the most massive (or 
the few most massive) cluster(s). Second, most of the star 
formation occurs near the end of the GMCs lifetime. We 
expect that this rapid burst of star formation is stopped by the 
disruption of the host GMC. 

Finally, the star formation rate is very different in GMCs 
than one would expect in models of star formation which es- 
sentially extend the Kennicutt relation (Equation valid 
on galactic scales, to the scale of individual GMCs, and as- 
suming that the star formation rate follows a Poisson process 

5 In Figure f3] we have not included a R c i = curve because it is nearly 
traced out by the S c j = 0.1 curve. 



over the entire cloud given by Equation (0 with random val- 
ues of T] that average to t) ~ 0.02, and where the cluster mass 
is determined by the random confluence of turbulent statistics. 

Rather, in the scenario proposed here, turbulence gener- 
ates clumps which subsequently grow by accretion. These 
clumps then experience runaway growth until stellar feedback 
ultimately unbinds the GMC, arresting the growth of all the 
clumps. This story highlights the importance of feedback on 
capping the ultimate star formation efficiency of roughly 10% 
in the Milky Way. 

To make observational contact, we calculate distribution of 
observed star formation rates in GMCs in our model. The 
star formation rate in a given GMC is driven primarily by the 
age of that GMC. Given a collection of GMCs with randomly 
distributed ages, the observed SFR distribution is simply the 
star formation rate as a function of GMC age convolved with 
the distribution of GMC ages. Very roughly, the likelihood 
of finding a GMC with a star formation rate in a given range 
is proportional to the length of time the most massive clump 
spends in the relevant mass range. 

In the upper panels of Figure [4] we show the histogram of 
star formation efficiencies per free-fall time 7790 for 779,0 > 
10 -4 . For (3^0 (middle and right panels), we show the cases 
for R c \ = 0.1 (blue solid line), 0.3 (green dotted line) and 0.5 
(red dashed line). These panels shows that a GMC spends 
most of its life at a very low 7790 ^ 0.02, but a significant 
portion of its life is spent at larger 7790 ~ 0.1. This is seen 
most clearly in the cumulative probability distribution in the 
lower panels of Figure |4] For instance, in the (3 = case, half 
of the lifetime of a cloud is spent at 7/90 < 0.015. However, 
10% of the lifetime of the cloud is spent at 779,0 > 0.1, with a 
maximum 779,0 =0.16. 

For more centrally concentrated GMCs, the fraction of 
GMC lifetime spent at large 779,0 become more pronounced. 
Indeed for, (3=1 and r = 0. 1 , the GMC spends half its life with 
?7ff,G > 0.08, and ~ 10% of its life is spent with 77^0 ~ 30%. 

The wide range in 7790 of the Bondi accretion model is con- 
sistent with the observed r ange in 7/9 q = 0.002 to 0.2 in Galac- 
tic GMCs dMurravH20iTh . Moreover, iMurravl (1201 lb found 
that GMCs with 779,0 ~ 0.2 are in the act of being disrupted. 
This is consistent with the notion that stellar feedback at high 
779,0 caps both the maximum star formation rate and the life- 
time of GMCs. 

We note that these observations are less consistent with 
models where global turbulent statistics set the SFR (Padoan 
1995; Krumholz & McKee 2005); in these models e ff = 77 = 
Vfffi ~ 0.017. Namely, for a turbulence-limited SFR, the val- 
ues of 77 can be expected to be roughly a factor of 3 above and 
below the average value of ij ss 0.02. However, the measure 
values of 7790 in our galaxy range up to w lOfff or higher. It 
is possible that this is merely a random conspiracy of turbu- 
lent statistics, but we would argue that accretion onto clumps 
gives a clearer picture of how such high SFRs are achieved. 

3.2. The effect of accretion on the mass spectrum of clumps 

The mass spectrum of both clumps and star clusters is gen- 
erally fit by a powerlaw distribution, 

dN d 



d In m 



/mo\ 
V m ) 



oc-\ 



(30) 



withal 1.5-2. e.g.. lKennicuttetalJ ( fT989l) . If, as we argue, 
clumps undergo a period of rapid growth, and that more mas- 
sive clumps grow more rapidly, will an initial power law dis- 
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FIG. 3. — Mass fraction, stellar fraction, and gas fraction of the clump as a function of time for fi = 1 (left), and 2 (right) for a GMC with Mqmc = 1O 6 M0. We 
show R c i = 0. 1 (black lines), 0.3 (blue lines), and 0.5 (red lines). For a clump growing near the center (r = 0. 1), note that the clump consumes the entire GMC is 
~ 1-5tqmc for P= 1 and 0.5tqmc for fi = 2. Clumps that are off center grow slower, but eve n slo west growing clump consumes the entire GMC in 3.5(1.5)tgmc 
for /3 = 1 (/3 = 2). For off-center clumps (r = 0.3 or larger), which, as noted below Equation {28) reduces to the fi = case, the clump experiences a long period 
of quiescience, followed by a period of rapid growth, but the period over which growth occurs is much shorter than in the fi = case. In these cases, we have not 
included a R c \ = curve because it is nearly traced out by the R c \ = 0. 1 curve. 





FIG. 4. — Upper panels: histogram of star formation rates over the lifetime of the GMC for fi = (left), 1 (center), and 2 (right). Values of rjgQ = 

r GMC^*/^GMC below 10~ 4 are discarded. Lower panels: Cumulative histogram of star formation rates rjn q over the lifetime of the GMC for fi = (left), 
1 (center), and 2 (right). The star formation rate is low for the majority of the lifetime of the GMC, but a sizeable fraction of a CMC's life is spent at rjg q > 0.1. 

That fraction increases with increasing central concentration, i.e., fi = 1 and 2. In both the upper and lower panels, we show the cases of R c [ = 0.1 (blue solid 
line), 0.3 (green dotted line) and 0.5 (red dashed line) in the middle (fi = 1) and right (J3 = 2) panels. 
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FIG. 5. — The cumulative clump mass distribution function N(> m). The 
red solid points, connected by a solid line, show the initial clump mass dis- 
tribution function, while the blue points (connected by a dashed line) show 
the clump mass distribution function at the point when the host GMC is dis- 
rupted. The initial slope is a = 1 .9, while the final slope is a = 1 .7. The initial 
total clump mass is 0.03Mqmc> the final total clump mass is 0.13Mgmc- 

tribution result in a final distribution of clumps (or star clus- 
ters) that is also a power law? 

Equation dT3l l shows how to map an initial clump mass into 
a final clump mass, given the age of the clump when accretion 
stops. We can map an initial clump mass function into a final 
clump mass function (when the host GMC is disrupted) by 
assuming a distribution of clump formation times. We work 
backwards from the time when the largest clump (or the com- 
bined effects of all the clumps) disrupts the GMC, when the 
largest clump has a given age. We then assume a Gaussian 
distribution of clump formation times with a mean equal to 
that of the largest cluster. To be definite, we assume that the 
largest cluster has a final mass of ~ 5% of that of the host 
GMC when the latter is disrupted. We also fix the number 
of clusters to be in the range 20-50, with total initial mass of 
O.OOIMgmc- We consider cases with 1.5 < a < 2. 

Figure [5] shows a typical result. The change in the power 
law slope a is Aa ~ 0.3, e.g., if a = -2.0 for the initial dis- 
tribution of clump masses, the final least squares fit (when 
the total stellar mass is ~ O.IMgmc ) is a = 1.7. More gen- 
erally, we find Aa 0.3 for a range in the inital slopes 
1.5 < a < 2.0. 

4. DISCUSSION 

4.1. Bondi Clumps in Numerical Simulations 

Recent simulations ( Vazque z-Semadeni et alJl2009l [2010) 
are consistent with our picture of star cluster formation. These 
authors studied the fragmentation of supersonic ally turbu- 
lent molecular gas into star-forming clumps (clouds in their 
nomenclature) and the subsequent evolution of these clumps. 
In their simulation, star-forming clumps occur when transonic 
converging flows in diffuse warm gas cause that gas to be 
rapidly transformed to cold dense material. This cold gas, 
which is in an initial pancake-like structure, collapses into fil- 
aments, and then to clumps. Their major result is that these 
star-forming clumps continue to accrete material from the 



background GMC at a high rate through the filaments. 

Acc retion via filaments i s also seen in smaller scale simula- 
tions dBaneriee et alj|2006h . In these simulations the collapse 
of a Bonner-Elbert sphere supported by supersonic turbulent 
motions is followed. The sheet-like and filamentary struc- 
tures that arise from the collisions of supersonic eddies lead 
to regions of very high density which collapse into protostel- 
lar clumps (disks and protostars). These authors found that 
the clumps continue to accrete material from the filaments at 
a very high rate, such that massive stars could form in a few 
thousand years. Indeed the accretion rate from filamentary 
accretion was a factor of 10 3 times larger than what would be 
expected from collapse of a singular isothermal sphere. This 
discrepancy arises from the assumption that the accretion pro- 
ceeds at the sound speed of the cold material; this sound speed 
is smaller than th e velocity needed to support an isothermal 
sphere of the type (Baneri ee et al.l20 06) simulated, since their 
clouds were supported by supersonic motions. 

In both sets of simulations just described, clumps are pro- 
duced in regions where transonic converging flows occur. 
Even in the absence of gravity from the newly formed clump, 
these convergent large scale motions will enhance the accre- 
tion rate of the clump relative to the rate that would be ex- 
pected from a similar mass clump placed at random in the 
simulation. We argue, however, that as the clump grows, the 
expanding gravitational reach of the clump will direct larger 
and larger amounts of gas to collapse onto the sheets and fila- 
ments, and thence onto the clump, yielding the accretion rates 
we have calculated. Testing this process by employing both 
gravitating and non-gravitating sink particles in large scale 
turbulence simulations is a subject of our ongoing work. 

4.2. The Effect of Vorticity 

So far we have ignored the possible role of angular momen- 
tum or vorticity in limiting the rate of accretion onto clumps. 
If the accreting gas carries a substantial amount of angular 
momentum, it may become rotationally supported at radii 
larger than the ~ pc scale at which we assume star forming 
clumps are born. If this were to happen, i.e., if the circulariza- 
tion radius r c ; rc of the accreting gas is larger than r c i, the rate 
of star formation would be suppressed relative to the rates we 
have calculated. 

We now argue that vorticity does not play a role in the for- 
mation of star clusters inside GMCs in local galaxies. 

Observational estimates of the angular momentum of 
GMCs start from mea surements of the veloc ity gradient Vv 
across the cloud, e.g., iRosolow skv et alj d2003l) . and assume 
£1 = |Vv|. As these authors note, strictly speaking this is an 
upper limit to f2, since | Vv| contains contributions from non- 
rotational motions. The measurements suggest that Sl(r) ~ 
const., i.e. the vorticity (V x v) is constant, or equivalently, 
that specific angular momentum j(r) of a parcel of gas at ra- 
dius r scales as j(r) = j(RcMc)(r / Rgmc) 2 ■ 

We characterize the rotation of the GMC by the parameter 



7 ; 



AgmcIVvI 
vgmc 



Then 



j(r) = 7/? gmc vgmc 



^GMC 



(31) 



(32) 



We assume that j does not change during the core accretion 
process, from its initial radius (the Bondi-Hoyle radius r%) 
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to the circularization radius r c ; rc , so that the circularization 
radius of a parcel of gas initially at radius r accreting onto a 
clump of mass M c \ is given by j(r) = V GM c \r c \ rc . We require 
that r c i rc < r c i, the cluster radius, leading to the critical value 



of the parameter 7 for a GMC with p(r) > 
CD) of 



r ' (c.f. Equation 



7c rit = ^ (1+W/2(3 -« 



Rcmc 



(33) 



If a GMC has 7 < jdu then gas will accrete onto a cluster 
of radius r c \ before the gas becomes rotationally supported, 
and the accretion can be described by Equation ( fTOb . Obser- 
vations of GMCs in M3 3 find typical values Rgmc ~ 23 pc 
dRosolowskv et al.l 120031) : scaling to /1 = 0.1 and (3 = 1, we 
find 7 cr i t = O.eeCO.l/^'^Crei/lpc) 1 / 2 . Typical values for 
the specific angular momentum of GMCs in M33 are j « 
22pckms _1 , corresponding to 7 w 0.17. Since 7 < 7 cl -i t it 
follows that clusters in M33 GMCs can grow to ~ 10% of the 
host GMC mass before vorticity begins to affect the accretion 
process. The situation is similar in M31, with Rgmc ~ 34 pc 
jRosolowskvll2007h . leading to 7^ = 0.55(0.1/^) 1/2 , while 
the observed value is 7 rj 0.28. Thus, we conclude that vor- 
ticity in GMCs is unlikely to limit the rate of accretion onto 
self-gravitating clumps, and hence the rate of star formation. 

5. CONCLUSIONS 

We have argued that star formation in turbulent GMCs, 
which is observed to occur in massive dense clumps, is 
in fact controlled by the properties of accretion onto those 
clumps. We start by assuming that the turbulence generates 
shocks, which form filaments and subsequently massive dense 
clumps. These turbulently generated clumps act as the initial 
conditions for the subsequent gravitationally dominated ac- 
cretion phase. We then used the observed properties of the 



internal gas motions in GMCs (Larson's law, vj(l) ~ l x l 2 ) to 
modify the Bondi argument for the accretion rate of a mas- 
sive dense clump embedded in that gas; for a clump located 
slightly away from the center of the GMC with a non-uniform 
density distribution, or for a uniform density GMC, that rate 
is given by Equation d29| >. Left to their own devices, these 
clumps can accrete all the gas in the GMC in only a few GMC 
free-fall times. 

We then showed that the star formation rate in an individual 
dense clump tends to track the mass accretion rate. The gas 
fraction in the cluster approaches a psuedo-fixed point, given 
by Equation ( f22l . 

The model predicts that the distribution of star formation 
rates in GMCs should peak at low values ijsfi < 0.01, but 
that there is a substantial tail to high star formation rates, 
with rjff c > 0.1, i.e., ten percent of the GMC is converted 
to stars in a single GMC free-fall time. This is in contrast 
to theories of turbulence regulated star formation, but agrees 
with recent measurem ents of % t Q on small spatial scales, e.g., 
Schru baet al.l (|2010), which show a broad distribution. A sec- 
ond prediction of the model is that the star cluster mass func- 
tion is flatter (has a smaller index a in Equation fl30l l) than the 
mass function of clumps, though it is still relatively close to 
the initial clump mass function. These aspects of our model, 
i.e., the broad distribution in star formation rates, the associ- 
ation of the largest star formation rates with the largest star 
clusters and GMC disruption, and a star cluster mass function 
with a power law slope of around a ~ 1 .5 - 2, are in broad 
agreement with observations. 

N.M. and P.C. are supported in part by NSERC of Canada. 
N.M. is supported in part by the Canada Research Chair pro- 
gram. This research has made use of NASA's Astrophysics 
Data System. 



REFERENCES 



Banerjee, R., Pudritz, R. E., & Anderson, D. W. 2006, MNRAS, 373, 1091 
Bigiel, F., Leroy, A., Walter, E, Brinks, E., de Blok, W. J. G., Madore, B., & 

Thomley, M. D. 2008, AJ, 136, 2846 
Bondi, H. 1952, MNRAS, 112, 195 

Bronfman, L., Cohen, R. S., Alvarez, H., May, J., & Thaddeus, P. 1988, 
ApJ, 324, 248 

Evans, II, N. J. 1991, in Astronomical Society of the Pacific Conference 

Series, Vol. 20, Frontiers of Stellar Evolution, ed. D. L. Lambert, 45-95 
Gao, Y, & Solomon, P. M. 2004, ApJ, 606, 271 

Grabelsky, D. A., Cohen, R. S., Bronfman, L., & Thaddeus, P. 1988, ApJ, 
331, 181 

Kennicutt, Jr., R. C. 1989, ApJ, 344, 685 
— . 1998, ApJ, 498, 541 

Kennicutt, Jr., R. C, Edgar, B. K., & Hodge, P. W. 1989, ApJ, 337, 761 
Klessen, R. S., Heitsch, F, & Mac Low, M.-M. 2000, ApJ, 535, 887 
Krumholz, M. R„ & McKee, C. F. 2005, ApJ, 630, 250 
Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 
Larson, R. B. 1981, MNRAS, 194, 809 

Leroy, A. K., Walter, F, Brinks, E., Bigiel, F, de Blok, W. J. G., Madore, B.. 

& Thornley, M. D. 2008, AJ, 136, 2782 
Li, P. S., Norman, M. L., Mac Low, M.-M., & Heitsch, F. 2004, ApJ, 605, 

800 

Mead, K. N, Kutner, M. L., & Evans, II, N. J. 1990, ApJ, 354, 492 
Molinari, S., et al. 2010, A&A, 518, L100+ 
Mooney, T. J., & Solomon, P. M. 1988, ApJ, 334, L51 
Mouschovias, T. C. 1976, ApJ, 207, 141 



Murray, N. 2009, ApJ, 691, 946 
— . 2011, ApJ, 729, 133 

Murray, N, Quataert, E., & Thompson, T. A. 2010, ApJ, 709, 191 
Murray, N, & Rahman, M. 2010, ApJ, 709, 424 
Padoan, P. 1995, MNRAS, 277, 377 
Rosolowsky, E. 2007, ApJ, 654, 240 

Rosolowsky, E., Engargiola, G, Plambeck, R., & Blitz, L. 2003, ApJ, 599, 
258 

Scheepmaker, R. A., Gieles, M., Haas, M. R., Bastian, N, & Larsen, S. S. 

2009, The Radii of Thousands of Star Clusters in M51 with HST/ACS, 
ed. Richtler, T. & Larsen, S., 103-+ 

Schruba, A., Leroy, A. K., Walter, F, Sandstrom, K., & Rosolowsky, E. 

2010, ApJ, 722, 1699 

Schruba, A., et al. 201 1, ArXiv e-prints 

Scoville, N. Z., Yun, M. S., Sanders, D. B., Clemens, D. P., & Waller, W. H. 

1987, ApJS, 63, 821 
Shu, F. H. 1983, ApJ, 273, 202 

Simon, R., Jackson, J. M., Clemens, D. P., Bania, T. M., & Heyer, M. H. 

2001, ApJ, 551,747 
Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 
Vazquez-Semadeni, E., Coh'n, P., Gomez, G. C, Ballesteros-Paredes, J., & 

Watson, A. W. 2010, ApJ, 715, 1302 
Vazquez-Semadeni, E., Gomez, G. C, Jappsen, A.-K., Ballesteros-Paredes, 

J., & Klessen, R. S. 2009, ApJ, 707, 1023 
Walcher, C. J., et al. 2005, ApJ, 618, 237 

Wu, J., Evans, II, N. J., Shirley, Y L., & Knez, C. 2010, ApJS, 188, 313 



2.5 



^ 2.0 

in 
"c 

3 

> 1.5 
03 

-M 

cr 

CD 

^ 0.5 



1 



0.0 



12 3 4 

SFR (arbitrary units) 



